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Abstract 

We study the possibility of determining the quark content of the scalar mesons a (980) and 
/° (980) through their hydrodynamical flow signature in relativistic heavy ion collisions. Utilizing 
the constituent quark scaling of hydrodynamic flow, we find that the tetraquark a (980) or / (980) 
mesons will have a of 0.38 at transverse momentum of 6 GeV/c in 20-60% central Au+Au 
collisions at t/snn = 200 GeV. The feasibility of measuring a°(980) — > tt t] and /°(980) — > tt + tt~ 
into the PHENIX and STAR detectors at the Relativistic Heavy Ion Collider (RHIC) is also 
discussed. Even though the mid-rapidity cross sections for these mesons at high-p-p are non- 
negligible, their broad mass range will make them difficult to detect in both p + p and Au+Au 
collisions. 
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I. INTRODUCTION 



Nearly all hadrons that have been measured can be described as a color singlet state 
of either three quarks or anti-quarks or a quark-anti-quark pair. However, the theory of 
Quantum Chromodynamics (QCD) allows more exotic structures to exist, like tetraquarks 
(qqqq) PQ, pentaquarks (qqqqq) |2j or multi-gluon states (glueballs) [3]. Several unexpected 
charm states have been discovered recently. For example, the X(3872) [I] and Z + (4330) [5], 
have not satisfactorily been explained as two-quark states but could be tetraquark states 

El El. 

Not only could exotic quark states exist at the higher masses, but several well-established 
low mass resonances are candidates for multiquark states. In Ref. PQ, Jaffe showed that the 
lowest mass scalar qqqq states form a nonet, several of which describe established particles. 
High-statistics measurements from radiative decays to the /°(980) |8j and a (980) [9] 
indicate that these are consistent with four-quark states. 

Assigning valence quark content for a particle is based on the determination of its quan- 
tum numbers. It would be extremely useful to have a more direct measurement of their 
valence quark structure. It appears that data from relativistic heavy ion collisions from the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory have yielded a 
dynamical observable that is sensitive to the number of valence quarks of a hadron. 

In collisions of gold nuclei at ^snn = 200 GeV, a quark-gluon plasma is formed and 
behaves as a strongly-interacting perfect fluid [TUl ITT] . In collisions with a non-zero impact 
parameter, the overlap region forms an almond-like shape. This is depicted in Fig.[TJ Because 
of thermodynamic pressure gradients in this region, particle production is not azimuthally 
symmetric. More particles are emitted along the reaction plane, the plane resulting from 
the beam and impact parameter axes, than perpendicular to it. This azimuthal asymmetry 
is characterized by the coefficients of a Fourier expansion of the invariant yield [12] 
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where pr, y, and <p are the particle's momentum transverse to the beam direction, rapidiy, 
and angle with respect to the reaction plane, respectively. The v n (pr, y) encode the strength 
of the azimuthal asymmetry. 
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FIG. 1. A beam view of a relativistic heavy ion collision. The impact parameter and the beam axis 
(into the page) form the reaction plane. The shaded, almond-shaped overlap region is the source 
of particle production, which is asymmetric in (p. 

When fluctuations in the initial energy and entropy density of the overlap region are 
small, symmetry requires only even Fourier coefficients to be non-zero. The dominant term 
in the series is v<±- Both the PHENIX and STAR experiments have measured a large number 
of identified baryon and meson t>2 as a function of pt [I3-16J. These data follow a universal 
trend in V2/n q vs. KEx/rig (see Fig. [2] and [3]). Here, n q is the number of valence quarks of 
the particle and KEt = mx — m = \jy\ + m 2 — m is the transverse kinetic energy of the 
particle. Since this data fits light-flavor (up, down, and strange) baryons and mesons for 
KEj-/n q < 1 GeV, a light-flavored particle's v<i is sensitive to the number of valence quarks 
of that particle. 

Such a result can be explained quantitatively by the recombination model of hadroniza- 
tion [T7j. When the quark-gluon plasma cools, hadronization must take place. The recom- 
bination model allows for quarks nearby in phase space to cluster and, with appropriate 
quantum numbers, form hadrons. The first model of recombination explained the measured 
enhancement of the p/ir + ratio in Au+Au collisions compared to p+p collisions. Because of 
the steeply falling momentum spectrum of quarks, the ability for three quarks of approxi- 
mately equal momentum combining to form a proton is more favorable than two combining 
to form a meson with the same final momentum as the proton. 
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In such a model, the quark v 2 will become the hadron v 2 based on the combination of 
the quarks that take place. The initial suggested scaling of v 2 was [T5] 

v 2 (p T ) ~ n q vl{p T /n q ) (2) 

However, detailed measurements of v 2 showed this scaling was incomplete and the KEt 
scaling was found to work better experimentally [15] . A more detailed recombination model, 
where one considers resonance-like scattering in the pre-hadronic state, which conserves 
energy and momentum in detail, is able to reproduce the observed KEt scaling [T9H2T] . 

With the discovery of constituent quark scaling of hadron v 2 at RHIC, we would like to 
search for evidence of particles with exotic quark structures. In this paper we present the 
study of the flow of light scalar mesons, / (980) and a (980), at RHIC. If the recombination 
model of hadronization is responsible for the translation of v 2 of quarks to v 2 of hadrons at 
RHIC, then these mesons will follow the observed scaling of v 2 and their number of valence 
quarks will then be known. A definitive measure of a tetraquark state would be an important 
result for QCD. 



II. FLOW OF SCALAR MESONS 

Using data it is possible to determine an expectation for the v 2 signature of the a (980) 
and /°(980). Fig. [2] shows the compilation of PHENIX data for v 2 of identified hadrons 
in 20-60% central Au+Au collisions at ^snn = 200 GeV [131 El- The v 2 /n q is plotted 
as a function of KEt = — m. Up to KEx/n q of approximately 1 GeV, all of the 
data follow the same trend. The PHENIX 7r° data were published in 10% centrality bins. 
They were combined for the figure using an eccentricity weighted average of the appropriate 
centrality bins. The eccentricity was determined by a Glauber Monte Carlo of the Au+Au 
collisions [22] • We fit a 4 th -order polynomial to the combined pion data: ir^ 1 for KEx/n q < 
1 GeV and 7r° for KE^jriq > 1 GeV. This is shown as the dashed line in Fig. [2J Fig. [3] is 
a compilation of identified hadron v 2 from STAR in minimum bias (0-80%) central Au+Au 
collisions at ^snn = 200 GeV [16]. A fit to the K s data, which exists over the broadest 
range in KEx/n q , is also shown as the dashed line. 

From here it is straight-forward to undo the scaling and plot the v 2 as a function of pt by 
knowing the mass and the number of valence quarks of the particle of interest. This is shown 
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KE^GeV) 

FIG. 2. (Color online) The combined PHENIX V2 data |X3|, I14j for identified tt^ (open circles), tt° 
(green closed circles), (red squares), p(p) (blue triangles), and <j> (magenta inverted triangles) 
in 20-60% central Au+Au collisions at ^snn = 200 GeV. The black line is a fit to the combined 
pion data. 
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FIG. 3. (Color online) The combined STAR vi data [16] for identified protons (green open circles), 
A (open squares), H (red crosses), ^ (magenta circles), and Kg (blue squares) in 0-80% Au+Au 
collisions at ^/snn = 200 GeV. The dashed line is a fit to the kaon data. 

in Fig. [4] for 20-60% collisions and Fig. [5] for 0-80% collisions for a particle of mass 980 MeV 
and consisting of either two (dashed line) or four quarks (dot-dashed line). There is a very 
large separation between these two curves. It is rougly consistent with the pocket formula 
Eq. [2j where the 4-quark V2 curve is roughly twice the 2-quark V2, especially at intermediate 
Pt where the recombination model dominates the particle production. In the mid-central 
collisions (20-60%), the 4-quark V2 at 6 GeV is 2.6 times larger than the 2-quark t>2 and 
reaches a value of 0.38, a v 2 that dwarfs all other t> 2 measurements from light hadrons. Even 
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FIG. 4. (Color online) The extrapolated v 2 of the a°(980) or /°(980) in 20-60% Au+Au collisions 
as a function of transverse momentum if they are either 2-quark states (dashed blue line) or 4-quark 
states (red dot-dashed line). For comparison the PHENIX pion v 2 is given as the solid (black) line. 
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FIG. 5. (Color online) The extrapolated v 2 of the a°(980) or /°(980) in minimum bias (0-80%) 
Au+Au collisions as a function of transverse momentum if they are either 2-quark states (dashed 
blue line) or 4-quark states (red dot-dashed line). For comparison the STAR v 2 is given as the 
solid (black) line. 

in minimum bias collisions the 4-quark v 2 reaches a value of 0.30 at 6 GeV and is 2.5 times 
higher than the 2-quark v 2 . 

We note that, in both the minimum bias and in mid-central collisions, the 4-quark v 2 
peaks between pt of 5-6 GeV or KEx/nq = 1-1.2 GeV. Therefore, the entire curve is within 
the measured KEt scaling of v 2 at RHIC. We expect the 2-quark v 2 to follow the trend of 
the 7T° or the Kg and decrease as a function of px, widening the gap between the expected 
v 2 signatures. 
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III. FEASIBILITY OF MEASURING SCALAR MESONS AT RHIC 



In order to establish that a measurement of the Vi can be made, we must get a sense of 
the signal strength. Important in v 2 measurements of resonances is the ability to measure 
the signal and the background Vi- In this section, we estimate the signal-to-background in 
p + p and Au+Au collisions to determine the feasibility of the measurement. 

In Au+Au collisions there will be at least three production sources: from radiative 
decays, "prompt" production from fragmentation, and production from recombination. We 
take the approach of estimating the first two components of these production mechanisms 
in p + p collisions. From this we can evaluate if the measurement in heavy ion collisions is 
feasible. 

Specifically, we are interested in the cross section of these mesons at high p?, i.e. at 
6 GeV where there is maximal difference between the 2-quark and 4-quark v^- We focus 
on the RHIC experiments where these mesons would be measured. The a (980) — > 7r°r/ 
decay would be possible in PHENIX where their calorimeter has measured both 7r° [23J and 
77 [21] inclusive production in p + p and heavy ion collisions. The /°(980) — > 7r + 7r~ would 
be a favorable decay for STAR, which has particle identification capabilities over a large 
kinematic range [25] . 

To estimate the cross section, we first study the radiative decays. PHENIX has mea- 
sured the cross section of and several other neutral mesons in p+p collisions at y/s = 200 
GeV The spectral shape for all mesons is well described with a Tsallis distribution [2T] . 

d 3 a 1 da (n - l)(n - 2) 

dp 3 2n dy (nT + m Q (n — l))(nT + m ) 
^nT + rnr 

K nT + m J K J 

where da/dy is the integrated cross section of the particle production at midrapidity. We 
modeled the radiative — > a (980)7 decay using the measured spectral shape from PHENIX: 
da/dy = 0.42 mb, n = 10, and T = 120 MeV. With the measured branching fraction from 
KLOE (7xl0~ 5 ), the da/dy for 6 GeV a°(980) production into the PHENIX pseudorapidity 
acceptance from decays is 9 nb. Requiring that the rj daughter decays to 2 photons, which 
is the most efficient method to measure 77 mesons in PHENIX, reduces the cross section to 
3 nb. If, instead, the would radiatively decay into an /°(980), based on the measured 
branching fraction from KLOE (10~ 4 ), da/dy for 6 GeV /°(980) into STAR would be 58 nb. 



If you require the /°(980) to decay into two pions, this reduces the cross section to 19 nb. 
These cross sections, therefore, represent the lower bounds on the production cross section 
for these scalar mesons into PHENIX and STAR, respectively. 

There should be an additional prompt production of a (980) and /°(980) mesons from 
fragmentation. This has not been measured at collider energies. We estimate the prompt 
cross section the following way. If the mesons are two-quark state, because their mass and 
quark content are similar to that of the if, their prompt production will be similar. PHENIX 
has measured the rf production cross section in p+p collisions [SH] . Its spectral shape is also 
given by the Tsallis distribution, Eq. [3] We use the PHENIX parametrization, da/dy = 0.7 
mb, n = 10, and T = 120 MeV, which yields a da/dy — 87 nb at 6 GeV. Requiring the r) 
from the a (980) decay into two photons reduces the cross section to 34 nb. The /°(980) 
decaying into two pions reduces the cross section to 29 nb. 

For a lower bound estimate of the prompt cross section, we use the following argument. 
If the mesons are 4-quark states, their production compared to 2-quark fragmentation would 
go approximately as the (p/vr + ) 2 . This is because this ratio encodes the difference in the 
string fragmenting into di-quark-anti-di-quark compared to quark-anti-quark pairs. In p + p 
collisions, this ratio is approximately 0.1. So we would expect that a lower estimate of the 
prompt cross section would be a factor of 100 less than the 2-quark production. In this case, 
the cross section for 6 GeV a°(980) into PHENIX and 6 GeV /°(980) into STAR is 0.34 nb 
and 0.29 nb, respectively. It is interesting to note that because the rf cross section at 6 GeV 
is the same order of magnitude as the cross section from decays, the lower bound values 
are significantly less than the radiative <fi production. Therefore, the decays would really 
represent an absolute lower bound for their production in p + p collisions. 

To find the signal-to-background ratio, we generated PYTHIA p+p events at y/s = 200 
GeV with Tune A [28J. For the a (980) background all 7c°-r] pairs, each decaying to two 
photons, were combined. Since this will be measured into the PHENIX acceptance, we also 
applied pt cuts consistent with the PHENIX capabilities [231 Elj- We required 7r° p T > 1 
GeV, rj pt > 2 GeV, and \y\ < 0.4 for both particles in the signal and background. These 
cuts reduced the signal cross section to 22 nb. The signal peak was parametrized from the 
KLOE data and normalized to have this cross section. The combined signal and background 
mass distribution is shown in Fig. [6j If the a (980) would have a cross section similar to 
the rf, then a peak is barely visible in the invariant mass distribution of 7r° — r] pairs. If one 
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FIG. 6. The 6 GeV pair pr tt ?? invariant mass distribution for background from PYTHIA (dashed) 
and the parametrized a (980) signal (solid) assuming the same cross section as the rj . In both the 
signal and the background, the ir° and rj decay to two photons and are required to be within the 
PHENIX y acceptance. The inset is the signal-to-background ratio over the full mass range of the 
a (980). 

integrates over the full mass range of 0.7-1.0 GeV, the signal-to-background is 10 -3 . The 
inset in Fig. [6] shows the signal-to-background as a function of invariant mass. It reaches a 
maximum of 3 x 10 -3 near 980 MeV. It follows that if the (f) decay is the dominant production 
mechanism, the signal-to-background will be reduced. With the pt cuts above, the signal 
from decay reduces to 2 nb over the whole mass range and the signal-to-background would 
be reduced to 5xl0~ 5 . Even though the cross section is small but measurable with RHIC 
luminosities, the fact that the resonance is spread over 300 MeV will make the measurement 
difficult. Additional cuts would be required to improve the signal-to-background. 

With the same PYTHIA background, we studied the /°(980) signal-to-background into 
the STAR acceptance. We assume the following cuts on the pion daughters of the /°(980) 
decay: pr > 1 GeV and \y\ <1. Even though the STAR experiment has particle identification 
capabilities below this, a higher-pr cut will naturally reduce the background. With these 
cuts the signal is reduced to 22 nb. The resulting signal and background 7r + 7r~ invariant 
mass distributions are shown in Fig. [7j Even though the /°(980) cross section into STAR 
is comparable to the a (980) into PHENIX, the 7r + 7r~ background is larger. An integrated 
signal-to-background is 4 x 10 -4 . The inset shows the signal-to-background over the entire 
mass range. 

With these upper estimates of the signal-to-background values in hand, we estimate 
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FIG. 7. The 6 GeV pair px tt + it~ invariant mass distribution for background from PYTHIA 
(dashed) and the parametrized J (980) signal (solid) assuming the same cross section as the rj . 
In both the signal and the background, the -k^ are required to be within the STAR y acceptance. 
The inset is the signal-to-background ratio over the full mass range of the /°(980). 

what they will be in heavy ion collisions. We assume the production of high-pr particles, 
4>, 7r°, rj and 7T , will scale as the number of nucleon-nucleon collisions, N co u, in the event. 
Therefore, we would naively estimate that the signal, which will scale as N co u, divided by 
the background, which scales as N^, would be reduced by a factor of N co u. In mid-central 
collisions (20-60%), a Glauber model yields an N co u of 200. This would be the reduction in 
the signal-to-background with no additional sources of these mesons in heavy ion collisions. 

There is an additional increase in the a (980) and /°(980) signal from recombination. We 
follow the estimation technique in Ref. [TTJ to obtain an order-of-magnitude estimate of the 
enhanced yield of tetraquarks due to recombiation. In that paper, the yield of mesons from 
recombination is related to products of the quark and anti-quark Wigner functions. For 
tetraquark states, the yield from recombination will be related to the products of two quark 
and two anti-quark Wigner functions. Assuming that the Wigner functions are exponential 
and independent of the quark or anti-quark internal quantum numbers, the a (980)/?/ will 
be a ratio of sums over the valence quark flavor, color, and spin. Using the notation of 
Ref. pi] this is 



where a, b, c and d represent the internal quantum numbers of valence quarks in the meson. 
This ratio is on the order of 10. The estimate assumes that the rj mesons will come predom- 



dN, 




(4) 



dN c 
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inantly from recombination. At 6 GeV, where we would like to apply this ratio, this might 
not seem to be a good assumption. However, as seen in Fig. [4] and [5j the tetraquark meson 
production mechanism will be dominated by recombination up to this point. We have also 
neglected the reduction of the rf production due to energy loss [21] . Including this estimate 
of the increase in yield due to recombination, the signal-to-background will decrease by a 
factor of 20 compared to p + p collisions. 

One could expect that there is additional signal from 7r° — 77 and n — n rescattering in the 
hadronic state, which would enhance the yield of these mesons and would mimic a 4-quark V2 
signal. We have checked this contribution using UrQMD [29] which simulates both a (980) 
and /°(980) resonance production in the hadronic state. We ran UrQMD in minimum bias 
mode and determined the absolute cross section for each resonance at 6 GeV. For this we 
used oauAu = 6.4 b [22]- For the a°(980) we determined a mid-rapidity cross section requiring 
the daughters to decay to photons and require the same kinematic cuts into PHENIX as 
we described above. We find a cross section of of 0.9 ph. This should be compared with 
the N co u- and recombination-enhanced cross sections fromx the the p + p production. From 
above the the and rf bound the cross section in p+p into the PHENIX between 2-22 nb. In 
Au+Au from N co u and recombination we expect the production cross section to be 4-44 ph. 
The lower bound is a factor of 4 above the hadronic-rescattering cross-section. Similarly, the 
/°(980) 6 GeV cross section into STAR using the same kinematic cuts as before was found 
to be 0.1 ph. The cross-section bounds on the /°(980) in p + p collisions is between 14-22 
nb. In Au+Au we would expect this cross section would be between 30-44 ph. Therefore, 
for /°(980) production, the hadronic-scattering component is negligible. For the a (980), it 
is a small component. 



IV. DISCUSSION 

From the previous section it is clear that the measurement of these resonances in heavy ion 
collisions will be difficult. Suitable background rejection must take place. An example would 
be placing higher pt cuts on the daughters in order to reduce the combinatorial background. 
However, increasing the 7r° cut to be greater than 2 GeV when considering a (980) decays in 
PHENIX, will only increase the signal-to-background by 30%. One might wish to search for 
the radiative <ft decays in heavy ion collisions. Requiring that the resonance combined with a 
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photon be in the mass range will certainly reduce the background significantly. However, 
the mesons from the radiative decays would not have the flow signature of interest. Because 
of the small Q of the decay, the a (980) or /°(980) would generally be aligned with the 
high-pr and exhibit the 2-quark v-i signature of the 0. Potentially, one could measure 
inclusive a°(980) and /°(980) and compare with those from decay since their vi signatures 
would be different. However, this would require a long running time to achieve the necessary 
statistics for such a measurement. 

This idea could be extended to other narrower tetraquark candidates, for example the 
X(3872) J/^7r + 7r _ . The advantage of measuring this state would be the narrow 3 MeV 
width [30]. The drawbacks for this measurement at RHIC is that the cross section is quite 
low considering the rate that was observed at Fermilab in inclusive production J3TJ 132] . 
Furthermore, because the charm quark is massive compared to light quarks, its t>2 is likely 
not the same as the light quarks [331 El]- The charm quark t> 2 has not been measured directly 
at RHIC, but only indirectly from the electrons from semi-leptonic decays [331 EH] • In that 
measurement there is a mixture of bottom and charm sources [371 138] . At increased RHIC 
luminosity, however, it may be possible to measure this resonance. Such a measurement 
might be more feasible at the Large Hadron Collider where the rate is higher. 

From the empirical evidence for constituent quark scaling, we expect a large V2 for the 
exotic scalar mesons a (980) and /°(980) in ^s^n = 200 GeV Au+Au collisions at RHIC 
if they are tetraquark states. Their v-i reaches a value of 0.38 in 20-60% collisions at 6 GeV 
and is a factor 2.6 higher than the 2-quark vi at the same pr- The mid-rapidity cross section 
of these mesons in p + p collisions will be larger than 19 nb for the /°(980) and 3 nb for the 
a (980), which results from their radiative decays. This is not a negligible cross section at 
RHIC. Because of the broad mass structure, the signal-to-background is at and below 0.001 
over the full 300 MeV mass range of the resonances. The signal-to-background will further be 
reduced in heavy ion collisions due to the increase in the underlying event multiplicity. Even 
with the enhanced production due to recombination, the signal-to-background will decrease 
by a factor of 20. Therefore, without taking steps to increase the signal-to-background, the 
measurement will be difficult. 
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